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Abstract 

This paper provides the connection between the Hankel transform and aerating trans- 
forms of a given integer sequence. Results obtained are used to establish a completely 
different Hankel transform evaluation of the series reversion of a certain rational func- 
tion Q{x) and shifted sequences, recently published in our paper [2]. For that purpose, 
we needed to evaluate the Hankel transforms of the sequences (a^C„ — l3Cn+i) and 
[a^Cn+i — ^^f^^, where C = {Cn)^^-^^ is the well-known sequence of Catalan num- 

bers. This generalizes the results of Cvetkovic, Rajkovic and Ivkovic [4]. Also, we need 
the evaluation of Hankel-like determinants whose entries are Catalan numbers Cn and 
which is based on the recent results of Krattenthaler [8] . The results obtained are general 
and can be applied to many other Hankel transform evaluations. 
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1 Introduction 

The Hankel transform of a given sequence a = {an)nmQ defined as the sequence h = (/in)^^^^ 
of Hankel determinants, i.e. 

hn = det {[ai+j]o<ij<n) , (n G No) (1) 

and denoted by /i = The term "Hankel transform" was first introduced by Layman [9] in 

2001. Despite that, many Hankel determinants evaluation were obtained much earlier, mostly 
due to their important combinatorial properties (see for example [3, 6, 13, 14]). 
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Papers [3, 4, 10] use a method based on orthogonal polynomials (or continued fractions) to 
provide a Hankel transform evaluation of different sequences. It is also used in our recently 
published paper [2] where we evaluated the Hankel transform of a series reversion of the function 
i+ax+i3x^ , as well as of the corresponding shifted sequences. 

In this paper, we also show another evaluation for the same sequences, which is based on 
the application of the falling a-binomial transform [12] and aerating transforms. The method 
described in this paper provides us with more general results regarding the connection between 
Hankel transforms and aerating transforms. 

For our purpose we need the Hankel transform evaluation of (a^C.„ — f3Cn+i)nQ^g 
{a'^Cn+i — /3C*.„+2)„gNQ, where C = (C*„)„gj!^^j is the well-known sequence of Catalan numbers. 
This generalizes results of Cvetkovic, Rajkovic and Ivkovic [4]. We also need the evaluation 
of Hankel-like determinants whose entries are Catalan numbers C„ and which is based on the 
recent results of Krattenthaler [8]. 

For our further discussion we need to recall the definition of the series reversion of a (gen- 
erating) function f{x) which satisfies /(O) =0 (see [1]). 

Definition 1.1. For a given (generating) function v = f{u) with the property /(O) = 0, the 
series reversion is the sequence (sn)„gNg such that 

u = f~^{v) = siv + S2v'^ H h H , 

where u = f~^{v) is the inverse function of v = f{u). Note that since /(O) = 0, there must 
hold So = f-\0) = 0. 



2 The series reversion of 



X 



\+ax+(3x^ 



We recall a few basic properties and expressions of the sequence (Mn)„gNQ obtained by reverting 

Q{x) 



1 + ax -\- 

This sequence is already investigated in [1,2]. The generating function U {x) satisfies Q{U {x)) 
X (Definition 1.1) and is given by 



- 1 - «x - a/1 - 2ax + (a^ - 4/3)a;2 

^^"^ = Wx • 

The general term of the sequence (un)^^^^ can be expressed in the following way (Proposition 
9 in [1]): 

E {'':Ac^o.-''-'p'- (3) 



Un = 

Note that the sequence (un)„gN(, generalizes the sequence (C„ + 5„o)„gNo a = 2 and /3 = 
1) and more generally the sequence {{Nn{,z) — 5no)/2^)„6No where Nn{z) is the n-th Narayana 
polynomial (for a = z + 1 and f3 = z). 

Consider the shifted sequences {u'^)neNo (""DneNo defined by m* = Un+i and m** = 
Un+2- Also denote by hn, and /i**, the Hankel transforms of the sequences m„, m* and m** 
respectively. Our previous paper [2] provides the evaluation of /i* , /i** and hn using the method 
based on orthogonal polynomials [4, 10]. The main results are the following theorems (Theorem 
4.3, Theorem 4.4, and Corollary 5.4 in [2]): 
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Theorem 2.1. [2] The Hankel transform of the sequence ('uJ^)„gp^Q is given by 

h: = ^("t'). (4) 

Theorem 2.2. [2] The Hankel transform of the sequence is given by 

h*: = ^/ Ua + v/«' - 4/3)"+' - (a - - 4/3)"+'] . (5) 

Theorem 2.3. [2] The Hankel transform of the sequence (Mn)„gNQ is given by 

hn = ^^'^ \(a- Ja^ -ApY - (a + Ja^ -ApY] . (6) 



Note that the sequence {u*^)nmo ^^duces to the sequence (C*^)„gpjQ of aerated Catalan num- 
bers (A126120), by choosing a = and (3 = 1 (yields directly from (2)). The sequence (C'")„gf^^ 
is defined by 

Cn/2, n is even 
0, n is odd 



c: 



a 



Also note that the Hankel transform of the aerated sequence (C!^)^^^^ is (l)„gpjQ (which is 
proven in the section 4), the same as in the case of the Catalan sequence (see for example [8]). 
That result raises the more general question about the Hankel transform of aerated sequences, 
which we deal with in the rest of the paper. 



3 The falling a-binomial transform 

The following transform is a generalization of the well-known binomial transform and was 
introduced by Spivey and Steil [12]. We will use it in further considerations. 

Definition 3.1. For a given sequence a = (an)„gN,)' ^^■^ falling a-binomial transformation 
b = B{a; a) is defined by 



n 

' n 



k=0 



Spivey and Steil [12] proved that the Hankel transform is invariant under the falling a- 
binomial transform for arbitrary a. In other words, the following lemma is valid. 

Lemma 3.1. [12] For an arbitrary sequence a = [a^njutma (md number a, itholds that'H{B{a]a)) 
H{a). 

The falling a-binomial transform can be written in the following matrix form 

h = B"a, B" = 



n,fcGNo 



where we treat the sequences a and h as corresponding column vectors (we also use this notation 
in the rest of the paper). We call the matrix B° the a-hinomial matrix. The following lemma 
shows the connection between the Hankel matrices 

and the matrix B". 

Lemma 3.2. If b = B{a\ a) then there holds 

= B"Ha(B")^. (7) 

Proof. Let us start from the general element bn+m of the matrix Hfe: 
Using the well-known identity 

\ t )=Z^\k)\t-k) 

we obtain 

n+m n / \ / \ m n / \ / \ 

t=0 k=0 \ / \ / 1=0 k=0 \ / V / 

m n ^ \ ^ \ m. n 

= E E ( 1) ■ ■ ( 7 ) = E E ■ • (B-).,., . 

1=0 k=0 ^ ^ ^ ^ «=0 k=0 

This completes the proof of the lemma. □ 

In the following sections, we give alternative proofs (to those given in [2]) of theorems 
regarding the Hankel transform of u*, u** and u„ (Theorem 2.1, Theorem 2.2 and Theorem 
2.3). 



4 The sequence u*^ 

We give an evaluation of the Hankel transform of n* just using transformations and known 
results concerning the Hankel transform of the Catalan numbers [8]. For this purpose, we 
define the following aerating transform. 

Definition 4.1. For a given sequence c = {cn)^^^^, we define its aerating transformation 
p = A{c) by 

{Cn/2, n is even 
0, n is odd 

In other words, if p = A{c) then p = (cq, 0, Ci, 0, C2, 0, C3, 0, . . .). 

Hence = A{C) where C = (C„)„gpj^ is the sequence of Catalan numbers. The following 
theorem shows the connection between the Hankel transform of a given sequence c and its 
aerated sequence p = A{c). 
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Theorem 4.1. Let g = 'H(p) and h = 'Hie) where p = A{c) is the aerated sequence of c. Then 
there holds 

det[pi+j]o<ij<n = det[cj+j]o<ij<[2.j ■ det[ci+j+i]Q^^j^^n^y 
In terms of Hankel transforms, this last equality can be written as 

9n = ^LtJ^[ilfij' 

where h* is the Hankel transform of the shifted sequence c* = (c„+i)^gj^y and h* = ^^(c*). 
Proof. By exchanging the rows and columns of the determinant det[pj+j]o<jj<n-i we obtain 

Co Ci 02 ■ ■ ■ 

ci C2 

Ci C2 C3 
'i+j\0<id<n = C2 Cg 
C2 Cs C4 



det [p; 



det 



B 



where A = [cj+j]o<jj<Lf J-i B = [cj+j_|_i]Q<^j<^n^j_^. Now the statement of the theorem 
follows immediately. □ 

It is known (see for example [7]) that "H ((C„)„gj^J = H ((C„+i)„gNo) = (l)neNo- Using 
Theorem 4.1 (for c„ = C„) we obtain the result V. ((Cn)neNo) = (l)neNo- 
We also need the following proposition: 

Proposition 4.2. Let c = (cn)^^^^ be an arbitrary sequence and h = 'Hie) its Hankel transform. 
Then H (('""c„)^gj^^) = {''^"'^"''^^^ hn) ^^^^ where r is an arbitrary number. 



Let On = P^Cn and let p = A{c). Recall that m* can be expressed as follows (directly from 



(3)): 



W„ — Un+l 



k=0 ^ 1=0 ^ ^ 



which implies that (M^)„gNg = B {p; a). 

Note that the sequence p = A{c) can be expressed as p = {(^^^'^C^)^^^^. Now using Lemma 
3.1 and Proposition 4.2 we obtain the result of Theorem 2.1: 



5 The Hankel transform of a linear combination of Cata- 
lan and shifted Catalan numbers 



Cvetkovic, Rajkovic and Ivkovic [4] considered the Hankel transform of the sequence (C„ + C„+i)^gj^^. 
In this section, we generalize their result, providing the Hankel transform evaluation of the se- 
quences {a'^Cn — /3Cn+i)ni^fq^ cind (tt^C^+i — f3Cn+2)neNo- uccd this result in the next 
section. 
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We use the method based on orthogonal polynomials, as used in [4, 10]. It assumes that a 
given sequence (a„)^gj^^ is a moment sequence with respect to some weight function (measure) 
wix), i.e. that there holds 

x"w{x)dx. 

If h = 'H(a) and /i„ 7^ for all n G Nq, then there exists a sequence of monic orthogonal 
polynomials (vr„(x))jjgj^^, which satisfies a three-term recurrence relation 

7r„+i(x) = (x - a„)7r„(x) - /3„7r„_i(x). (8) 

The Hankel transform hn can be evaluated using the following Heilermann formula (see for 
example [7]): 

hn = a^^'(3^(3r'---(3Li(3n. (9) 

In order to establish closed- form expression for coefficients an and Pn, the following transfor- 
mation lemmas can be useful: 

Lemma 5.1. [2] Letw{x) andw{x) he weight functions and denote by (vr„(x))^gj^^ and (vrri(a;))^gj^^ 
the corresponding orthogonal polynomials. Also denote by (a„)„gj^^ , (/3n)„gNo ^^'^ ('^")neNo ' ('^") 

the three-term relation coefficients corresponding to w{x) and wlx) respectively. The following 
transformation formulas are valid: 

(1) If w{x) = Cw{x) where C > then we have an = oin for n G Nq and (3^ = Cf3Q, f3n = /3n 
for n G N. Additionally there holds 7r„(x) = '/r„(x) for all n G Nq. 

(2) If w(x) = w{ax + h) where a, 6 G M and a ^ there holds an = for n G No and 
^0 ~ fa\ ^'^^ ~ f^'^ n G N. Additionally there holds vr„(x) = -^7rn{ax + b). 

Lemma 5.2. (Linear multiplier transformation) [5] Consider the same notation as in 
Lemma 5.1. Let the sequence (r„)^gjs^^ be defined by 

ro = c-ao, rn = c-an — (n G No). (10) 

rn-i 

If w{x) = {x — c)w{x) where c < inf supp(u'), then there holds 

~ f ~ T 

/3o = / w{x) dx, Pn = Pn—^, [n G N), 
an = ttn+i + rn+1 - r„, (n G No). 



n£No 



Now we prove the main theorem of this section. The proof is based on the sequential 
application of the previous two lemmas for a known weight function, i.e. a weight function 
whose coefficients a„, and /3„, are known. 

Theorem 5.3. The Hankel transforms of the sequences {a^Cn — PCn+i)n(zna ^'^^ {oi^Cn+i — PCn+2)nmo 
can be evaluated as follows 



det[a Cj+j — /3Ci+j+i]o<jj<n 



1 



detfa^Ci+j+i — pC, 



22n+3^C^2 -4/3 L 

1 



{a + ^a^ - 4/3)2"+3 -{a- ^a^ - 4/3)2"+3 



i+j+2j0<i,j<n 
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Proof. We use again the method based on orthogonal polynomials. It is well-known (see for 
example [4]) that the Catalan numbers have the following moment representation 



1 r'. 14, 



This directly implies that the sequence (a^C„ — (3Cn+i)n(zNo moment sequence of the 

weight function: 

1 /4 
Zir \ X 

In order to determine the three-term recurrence relation coefficients corresponding to we 
start from the weight function of the monic Chebyshev polynomials of the fourth kind: 



The corresponding coefficients an^ and (5n^ are 

= -1/2, af = 0, n>l, C = 7r, = 1/4, n>l. 

Now we define a new weight function w^^\x) by 

and use part (2) of Lemma 5.1 with a = 1/2 and b = —1. Hence we obtain 

4') = 1, «i') = 2, n>l, tf^=27r, /^i^^ = 1, n > 1. 
The next transformation is 

w^^\x) = -^-w^^^x). 

ZTT 

From part (1) of Lemma 5.1 we see that 

af = l, ai') = 2, n>l, = -/?, Pj^^ = I, n > 1. 

The final transformation is a linear multiplier transformation 



,{2)/ 



a 

w(x] = \ X — ^ ' ' ""^'''t^'- 



According to Lemma 5.1 (c = / (i) we have to consider the following temporary sequence 
and the coefficients /3„ are obtained by 

rn-i 

The Heilermann formula now yields 

hn+l 
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Replacing the last expression into (12) we obtain the following linear difference equation 

- («' - (3)hn-i + (3%,.2 = 0, {n> 2) (13) 

where the initial values are given by /iq = — /? and hi = a'^ — 3q;^/3 + By solving (13), 
we directly obtain the first statement of lemma. 

To prove the second statement, let us observe that the weight function of the sequence 

(a'^Cn+l - (3Cn+2)nsNo ^^^^^ 



2 



w(x) = — x(a — 0x)\ 1 = — ■ I — X 

The initial weight function now is the weight function of Chebyshev polynomials of the second 
kind: 

u;(o)(x) = Vl a;G[-l,l]. 
As in the previous case, by applying the following sequence of transformations 

«;(i)(x) = ( , w^^\x) = ■ w^^\x), w{x) = fa; - ^^^\^) 



we prove that the sequence \hn) satisfies the same linear difference equation 

K - («' - P)K^i + = 0, {n> 2) (14) 

but with different initial values = — 2/3 and hi = a'^ — 3a'^f3 + 3/3^. By solving (14) we 
obtain the second statement of lemma. □ 



6 The sequence u 



n 



We can also express the sequence «** = ^^+2 as the falling a-binomial transformation of a 
certain sequence, whose Hankel determinant will be evaluated. First we need to define the 
generalization of the aerating transform A{p). 

Definition 6.1. For a given sequence c = {cn)neno' '"'^ define its a-aerating transformation 
b = A{c;a) by an = apn + Pn+i, where p = A{c). In other words, if a = A{c;a) then 
a = {aco, ci, aci, C2, ac2, cs, acs, . . .). 

Let a = A{c;a), i.e. the a-aerating transform of the sequence c„ = /3"'Cn- The sequence 
a = (an)„gNo '^^^ expressed as follows 

a(3''Ck, n = 2k 
^ P^Ck, n = 2k-l 

According to (3) we have 

rn + li 



(15) 



k=0 ^ ^ 

fc=0 ^ ^ A:=0 ^ ^ 

n 



1=0 



^ n—l „ 

a ai- 



Hence u** = B{a; a) and from Lemma 3.1 we conclude that ^{u**) = 'H{a). To evaluate 'H{a), 
we need the following theorem. 

From now on, we denote by [Aj^xm a matrix formed by first m rows and columns of the 
(infinite) matrix A. Also we label rows and columns of matrices starting from (i.e. 0, 1,2,...). 

Theorem 6.1. Let g = 'H{a) and a = A{c] a). Then the following holds 
9n = det[aj_|.j]o<ij<n 



det [a^c. 



i+j Q+j+l 



|o<ij<fc-i ■ det[Q+j+i]o<jj<fc-i, n — 2k — 1 



(16) 



a ■ detfa^Cj+j+i — Cj+j4.2]o<i,j<fc-i ■ det[cj+j+i]o<ij<fc, n — 2k 
Proof. The determinant det[aj+j]o<ij<n that we wish to evaluate, has the form 



det \a 



0<i,j<n 



det 



aco c\ aci C2 

Ci aci C2 ac2 

aci C2 ac2 C3 

C2 ac2 C3 ac3 



J (n+l)x(n+l) 



We distinguish two cases depending on the parity of n. 

Case 1. = 2A; — 1 is odd. We multiply column 2j + 1 by and subtract from the column 
2j, for every j = 0, 1, . . . , A; — 2. That leads to the following determinant 



det [ai+j,]o<i,j<n 



det 



acQ — a ^Ci C\ aci — a ^C2 C2 





- c 




aci 

C2 

ac2 



aci — a ^C2 C2 ac2 — a ^03 C3 





- c 




ac2 
C3 



(n+l)x(n+l) 



By exchanging rows and columns we get the block diagonal form: 



det \a 



i+j \0<i,j<n 



det 



A * 
B 



det A ■ detB, 



where star (*) denotes the appropriate k x k matrix which does not have an influence in 
determinant computation. Matrices A and B are given by 



[aci+j — a Q+j+iJo<jj<fc-i, 



B 



aa 



'j+j+lj0<i,j<fe-l- 



By taking a from each column of matrix B and putting to the corresponding column of matrix 
A, we get the first case of the expression (16). 

Case 2. n = 2k is even. Multiplying column 2j by and subtracting from the column 2j — 1 
(for every j = 1,2, . . . , k) yields the determinant 



det [a 



det 



a Co 







aci 







ac2 


Cl 


aci 


- a~^C2 


C2 


ac2 


- a"^C3 


C3 


aci 







ac2 







ac3 


C2 


ac2 


- a^^C3 


C3 


ac3 


— a~^C4 


ac3 


ac2 














ac4 



J (n+l)x(n+l) 
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Again, by exchanging rows and columns we get the block diagonal form: 



det[ai+j]o<ij<n = det 

where 



A' * ■ 
B' 



det A' • det B', 



A' — [«Cj+j-|.i]o<jj<fc, B' — [acj_|_j_|_i — a ^Ci^j^2\o<i,j<k-i- 

By taking a from each column of the matrix B' and putting the first k entries to the corre- 
sponding columns of the matrix A', we obtain the second part of (16). □ 

According to the previous theorem, Proposition 4.2 and the fact that "H {{Cn)nmo) ~ 
(-'-)nGNo' '^^^ n = 2k — 1, it holds that 

det[aj_,.j]o<jj<„ = det[a^/3*"'"-'Cj+j — (3^^^~^^Ci+j+i]o<ij<k-i ■ det [/?*"'"■'"'" "'"Cj_|_j_|_i]o<jj<fc-i 
= /^(a) det[a'^Ci^j — /3Ci+j+i]o<ij<fc-i ■ det[Cj+j+i]o<jj<fc-i 
= /^(a) det[a'^Ci+j - /3Ci+j+i]o<ij<fc-i. 
Similarly, for even n = 2k, we have 

det[ai+j]o<i,j<n = adet[a'^f3'+^+^Ci+j+i - f3'+^+^Ci+j+2]o<i,j<k-i ■ det[/3'+^+^Ci+j+i]o<i,j<k 
= a (3" ^^"^^ det[a^Ci+j+i - f3Ci+j+2]o<i,j<k~i- 
Now expression (5) (i.e. Theorem 2.2) is directly obtained using Theorem 5.3. 

7 Hankel-like determinants based on Catalan numbers 

Before we proceed to the reevaluation of the Hankel transform of the sequence ('Un)„gNo' 
need to prove two lemmas concerning determinants which are generalizations of the Hankel 
determinants. Our main tool is the following theorem proven by Krattenthaller in [8] (Theorem 
3): 



Theorem 7.1. [8] Let n be a positive integer and ao, ai, . . . , non-negative integers. Then 

{i + k)\{2a,)\ 
{2i)\ai\{ai + k)\' 



k—1 ^ A^^I^Sq!'^! 

det[c„,+,]o<,,<,_i = n K-«^)n J^u.ir^.i^Mi - (17) 



0<i<j<k-l i=0 



We use the notation x{P) = 1 if -P is true and xiP) = otherwise. Also, we assume that 
the sequence c — (c„)^gj^^ is defined by c„ — P^Cn- 

Lemma 7.2. For every / = 0, 1, . . . , /c — 1 we have 



A^ilr 1 _ okHk-ifl + k + l\ 

aet[Ci+j+xQ>i)+ijo<ij<fe-i - p \ 21 + 1 )' 
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Proof. Denote by = i + x(i >/) + !. We need to evaluate the determinant 

det[ca^_i_j]o<ij<fc_i = detlP"'-^^ Cai^j]o<i,j<k-i- 

By taking from row i [i = 0,1, . . . , k — 1) and then from column j [j = 0,1, . . . , k — 1) 
we get the total power of f3 equal to 

k-l k-l 

+ =k^ + k-l. 

Hence, our determinant reduces to 

det[Ca,+j]o<i,j<k-l = Z?'^ '^'^ ^ det[Ca,+j]o<i,j<k-l- (18) 

According to the Theorem 7.1, we need to compute the following products 

Pi= n ^2=n 



0<j<j<fc-l i=0 ^ ' « V ' ' ^ 



By direct evaluation it can be shown that 



n 



The first product was evaluated considering three different cases {i<j<l,l<i<j, 
i < I < j), while for the second we only needed to distinguish between i < I and i > I. Now 
using (19) and Theorem 7.1 we obtain 



det[Ca,+j]o<i,j<k-i = Pi - P2= { 21 + 1 / 



Now the statement of the lemma follows directly from the previous equation and (18). □ 
Lemma 7.3. For every I = 0,1, . . . ,k — 1 we have 



det[cj+j+^(j>;)]o<ij<fc-i — (3' 



k^~i 



l + k 
21 



Proof. We use again the same procedure. Now we denote = « + x(^ ^ obtain 

det[Ca,+j]o<i,j<fc-l = det[/3°''~^^Cai+j]o<i,j<k-l = ^ det[Cai+j]o<i,j<k-l- 

The product Pi has the same value as in the previous case, while P2 is equal to 

(/ + 1)! il + k)\ 

j=0 

Again, by replacing det[Cai+j]o<i,j<k-i = Pi ■ P2 (Theorem 7.1) we obtain the statement of the 
lemma. □ 
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8 The sequence Un 

In the section 6 we proved that u** = B{a]a), where a is the a-aerating transform of the 
sequence Cn = (3"'Cn {a = A{c; a)), i.e. (equation (15)): 

a(3''Ck, n = 2k 
f3''Ck, n = 2k-l' 

According to Lemma 3.2, we have H^j.. = B"H(i(B")'^. We have already proved that u* = 
B (p; a) (section 4) where p = A{c), i.e. 



Pn 



l3^Ck, n = 2k 



0, 



n = 2k-l 



This can be written in matrix notation as u* = B"p. Now we have that the following matrix 
equality holds: 



'1 




'0 


p 




'1 








_p 


Ha 









p^iB'^f 
B> B°H„(B"^^ 





'0 











H, 



(20) 



Hence, the determinant of the (n + 1) x (n + 1) principal minor of H^, formed by the rows and 
columns with indices l,2,...,n + l, is equal to the same minor of the matrix 



H' 



p^ 

P Ha 



That minor is exactly hn, i.e. ?7,-th member of the Hankel transform h = 'H(u). In other words, 
we have to compute 



hn = det[H„](„+i)x(n+i) = det 





P Ha 



det 



•^J {n+l)x{n+l) 






Co 





Cl 





Co 


aco 


Cl 




C2 





Cl 


aci 


C2 




Cl 




C2 




C3 





C2 


ac2 


C3 


ac3 



J (n+l)x{n+l) 



We distinguish two cases depending on the parity of n. 



Case 1. 77, = 2 A; — 1 is odd. Multiplying the column 2j by a and subtracting from column 
2j + l (j=0,l,...,/c-l) yields 



hn = det 






Co 







Cl 





Co 





Cl 







C2 





Cl 


aci 


C2 


— a^ci 


ac2 


Cl 





C2 







ca 





C2 


0!C2 


C3 


- «^C2 





J {n+l)x{n+l) 

By exchanging rows and columns in the previous determinant we obtain 



;-l)*^det 



A * 
B 



-l)Met A ■ detB 



(21) 
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where the matrices A and B are equal to 

Co Ci 



A = det 



Ci C2 — a Ci 



Cfc-l 
Ck-2 - oi^Ck-l 



B= k 



'i+jj0<j,i<fc-l- 



Ck-l Ck-2 - a Ck-1 C2fc_2 - a C2fc-3_ 

Since c„ = (3"'Cn, using Proposition 4.2 we obtain 

det B = det[(3'+^a+,]o<ij<k-i = (3'^'-'^ det[a+,]o<.,,<fc_i = f^'^'-'l (22) 

By adding column j to column j + 1 of the matrix A (j = 0, 1 , . . . , /c — 2) we obtain the following 
determinant 



det A = det 



Co a^co + ci a^co + a^ci + C2 



Cl C2 
C2 C3 



C3 
C4 



- fcxfc 



Expanding over the first row yields 



k-l 



det A = ^(-1)' I ^a^'^'ci-h j (^'si[ci+j+xij>i)+i\o<i,j<k-2- 

1=0 \h=0 J 

Using Lemma 7.2 together with the expressions (21), (22) and (23), we obtain 
h2k-i = (3^'-'^^''-'^J2(-^^'^' 



1=0 



^h=0 



(23) 



(24) 



Case 2. n = 2k is even. Subtracting column 2j — 1 from column 2j (j = 0, 1, . . . , /c — 1) we 
obtain 

Co — aco Cl 



h„ = det 



Cq aco Cl — a^Co aci C2 — a^Ci 



Cl 
Cl aci 

C2 





- 





Cl aci C2 — a^Ci ac2 C3 — a^C2 



C2 

ac2 
C3 





- 




J {n+l)x(n+l) 

By exchanging rows and columns in the previous determinant we obtain 



hn = (-l)*^det 



A * 
B 



-D'^det A ■ detB 



where the matrices A and B are equal to 



A = det 







— aco 



Cq Cl — a Co 



-acfc_i 
Ck - a^Ck-i 



B — [Ci+j+l\o<i,j<k-l- 



Ck-l Ck - a^Cfc_i C2fe-i - a C2fc-2 

Since c„ = P^Cn, using Proposition 4.2 we obtain 

detB = det[/3*+^^+ia+,-+i]o<i,,<fe_i = (3''" detfC, 



_ rik^ 

+iJo<j,i<fc-i = P ■ 



(25) 



(26) 
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Again, by multiplying column j by a and adding to column j + 1 (j = 0, 1, . . . , — 1), we 
obtain the following determinant 



det A = det 






—aco 


a^co — aci 


— a^co — a^ci — ac2 


Co 


Cl 


C2 


cz 


Cl 


C2 


C3 


C4 


C2 


C3 


C4 


C5 



(fe+l)x(fe+l) 



which can be expanded by the first row in the following way 

k //-I \ 

detA = ^(-1)' I ^a^^'+^Q-i-h j det[ci+j+^(j>/)]o<jj<fc-i. 

1=1 \h=0 J 

Now using Lemma 7.3 and expressions (25), (26) and (27) we obtain 

k / l-l 



(27) 



h^k = ^(-1)'=+'-^ ^a^h+i^fc-i-/.^^ 



z=i 



,/i=0 



l + k 
21 



(2J 



Proof of Theorem 2.3. We can rewrite expressions (24) and (28) as follows (we exchanged 
the order of summation): 



k-i 



fc-i 



h=0 l=h 
h,, = ^M2fc-1) J2 (-1)^+'-!^ 



l + k 
21 + 1 

l + k 



(29) 



h=0 



l=h+l 



21 



Now let Zn = /3 ^'^^hn and in the second equation of (29) decrease the bounds for / by 1. 
Expressions for Z2k and Z2k-i are 



Z2k- 



h=0 l=h ^ 

k = Y.^'^^"(^'~"''Y.^-^)'^'^^A 

h=0 l=h ^ 



(30) 



l + k + l 
21 + 2 



By direct verification we conclude that Zn satisfies the three-term linear difference equation 

Zn+2 - OlZn+1 + (3Zn = 

for all n G No, which directly implies the expression (6): 

hn = P^^^Zn = ^^'^ 

2V«'-4/3 
This completes the proof of Theorem 2.3 



)Tl / \ Tl' 
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